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Coherent control of a quantum system is limited both by the decoherence due to environment and 
the quantum nature of the control agent. The high fidelity of control demanded by fault tolerant 
quantum computation and the intrinsic interest in nonclassical effects from the interplay between 
control and dissipation are motivations for a detailed study of the interaction dynamics between the 
quantum system and the macroscopic environment and control agent. We present a detailed time 
evolution study of a two-level system interacting with a laser pulse and the electromagnetic vacuum 
in the multimode Jaynes-Cummings model. A diagrammatic formalism allows easy identification of 
coherent dynamics and relaxation of the two-level system. We demonstrate a computational method 
of dynamics with precise error bounds for fast operations versus slow decoherence, spanning the 
Markovian and non-Markovian regimes. Comparison against an exact model solution of our results 
with existing approximations of the master equation shows the lack of accuracy in the latter. 
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I. INTRODUCTION 

An open system, i.e., a small quantum object in the 
presence of a macroscopic environment, presents a fun- 
damental problem in quantum mechanics and its appli- 
cations. We wish to address here the problem of the envi- 
ronment with the dual function of decoherence and con- 
trol of the quantum object. For practical purposes, both 
experiments on coherent processes and quantum technol- 
ogy require a small parameter to / Ti in the time scale of 
the control duration to being much smaller than the deco- 
herence time Ti . A paradigmatic system for this problem 
is the interaction of a two-level system (TLS) with the 
quantized electromagnetic fields^ In terms of the TLS- 
photon interaction strength g in units of frequency, the 
controlled TLS process (the Rabi rotation) is a strong 
coupling process with g|a|to ~ 0(1) using a coherent 
photon state \a) with a large mean number \a\ 2 of pho- 
tons while the long decoherence time is a weak coupling 
process t^/Ti <C 1 within the control time. In particular, 
the high fidelity of the operation to an error threshold 
between 10~ 3 and 10" 4 in fidelity demanded by fault 
tolerant quantum computation^ sets the bar for high 
accuracy in the theory of the open system. While the 
TLS open system problem has been much studied, we 
posit that the specific additive problem of the control 
and decoherence processes remains. 

The decoherence problem of a TLS in a spin bath and 
its suppression under a classical control have been exten- 
sively studied; 4 -^ but the noise due to a quantized optical 
control was not taken into account. Barnes and War- 
ren demonstrated the decoherence induced by the back 
action from the TLS to the electromagnetic control^ 
For a Markovian system, this problem can be solved 
using the optical Bloch equation^. However, in many 
non-Markovian systems with a structured environment, 
say nanocavities^ and photonic band gap materials^, 
the optical Bloch analysis no longer applies. Because 



of the large Hilbert space of the multimode TLS-ficld 
Hamiltonian, an exact diagonalization of the problem is 
impractical^ Therefore, a quantum theory suitable for 
a non-Markovian open system with high accuracy set for 
our problem is needed. 

In this paper, we develop a formalism to solve the 
multimode Jaynes-Cummings (JC) model under a coher- 
ent light pulse with arbitrary pulse shape in the limit 
to <C Ti, relevant to the problem posed. The exact non- 
unitary evolution of the TLS is expressed in terms a time 
evolution forward and a reverse evolution backward of 
emitted and absorbed photons conditioned on the TLS 
down and up transitions. The photon dynamics are eval- 
uated by the field theoretic perturbation series. The di- 
agrammatic structure allows an explicit identification of 
the photons' role in coherent dynamics and dissipation 
process of the TLS. All perturbation terms in a coher- 
ent dynamics segment are summed and the dissipativc 
processes are expanded in powers of a small parameter 
(t-o/Ti) 1 where 7 depends on the photon correlation time 
t c and ranges between 1 (the Markovian limit) and 2 (the 
non-Markovian limit). The process gives a practical nu- 
merical computation procedure for a given error limit in 
powers of the small parameter. Wc illustrate the compu- 
tation for a short-time quantum operation of the TLS by 
including all the relevant quantum processes within the 
first order error bound. We identify the precise origins 
of the interference effects between the control and the 
dissipation in the evolution processes. 

The master equation (ME) approach of treatment 
of the open quantum system has been enormously 
important for quantum opticsii and for quantum 
information^ We have made a comparison of our the- 
ory with several prominent approximations extant in the 
ME approach 1 ^ for a model problem involving the inter- 
ference effect between control and decay which has an 
exact solution and a semiclassical one. The compari- 
son results show that in the non-Markovian regime, the 
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driving forces for the coherent control and the dissipa- 
tion are not additive and that the ME approximations 
are closer to the semiclassical results than the quantum 
results, not meeting the stringent accuracy requirement. 
Perhaps, the comparison results would stimulate effort to 
refine the ME approximations. 

The outline of the rest of the paper is as follows. In 
Section|lIl we develop a field-theoretic solution for the JC 
model in the presence of the multimode photonic field, 
elucidating the coherent and dissipative components from 
the diagrammatic structure. In Sec. lIIIl a detailed analy- 
sis of the effect of vacuum decoherence on the control pre- 
cision of the TLS is given in the non-Markovian regime. 
Sec IIVI compares the diagrammatic solution with the ME 
approximations. Sec. [V] summarizes and the appendices 
add some technical details. 



II. TLS-PHOTONS INTERACTION DYNAMICS 

By the field theoretic techniques, we express the evo- 
lution operator of the whole system (the TLS and the 
photons) as an infinite perturbation series in terms of 
the spin and photon operators. Then, we evaluate the 
matrix elements of the spin exactly, resulting in a dia- 
grammatic series of the photon operators only. By the 
Wick's theorem on the photon operators, we build a per- 
turbativc solution to the non-equilibrium problem of the 
dynamics of the laser photons and the TLS in the bath of 
the electromagnetic vacuum. We find a controlling small 
parameter, (io/Ti) 7 , defined in the Introduction and de- 
tailed below, for the perturbation series. The approach of 
removing the spin operators first stands in contrast to the 
standard ME approach^ which traces out the photonic 
environment first and then solves the equation of motion 
of the TLS, and which lacks error bounds for most of its 
approximations. 



A. TLS transformation by coherent photon state 

We start with the canonical multimode JC 
Hamiltonianjiii^ 



H = H + V, 



where H = ]-uj <t z + ^ 



w k a\a k , 



g k \a+a k + &-a k 



(1) 



(2) 



Hq contains the bare Hamiltonian of the TLS of energy 
splitting ujq with the Pauli operator a z and the photons 
of energy oj k with creation operator a k . V is the TLS- 
photon interaction with the coupling constant g k , pre- 
sented in the rotating wave approximation (RWA) which 
is justified in Appendix [XJ While the single-mode JC 
may be used to treat the laser pulse by making the cou- 
pling time-dependent g(t), the multimode extension^ is 



more suited for our purpose of investigating the joint 
quantum effects of the light control of TLS and the dis- 
sipation. 

The composite system of the TLS and photons is 
given by a product initial wavefunction, |^(0)) = 
E s c sl s )] \ a )> wncrc s = ± denotes the two states of 
the TLS and a = (a kl ,a k2 , . . .) denotes the multimode 
coherent state. The reduced density matrix P Sft3 '(t) is 
expressed in terms of the transformation matrix condi- 
tioned on the initial photon state: 

P s f ,s' f (t) = 2J c sCs'Ps f , S ' f ; S , S '(t, a) (3) 
s,s f 

P.„S f = (cc\(s'\U\t)\s' f )(Sf\U(t)\s)\ a ) 

X e i ( s '/ 1 - s / 1 )"ot/2_ (4) 

In the interaction picture, the time evolution operator 
and the interaction are: 
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and A k = loq — uj k is the detuning of the k mode. This 
form of transformation of the reduced density matrix vis- 
ibly retains the quantum nature of the evolution of the 
composite system and is easily reduced to the problem 
of the expectation value of the electromagnetic field op- 
erators for the initial photon state. The initial product 
state may be generalized to any composite state |\P(0)) = 

J2 S J D«c s , Q \s) |a>, where / Da = / . . . 

In the perturbation series of the evolution operator, the 
TLS state is flipped up or down by a series of interaction, 
leaving only the corresponding photon operators, 
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B. TLS evolution in terms of photon processes 

The task is reduced to evaluate the transformation ma- 
trix in Eq. (jj) in terms of the series expansion of the uni- 
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tary operator and its inverse in Eq. (JTJ) or (JHJ that consist 
of photonic components in Eqs. The formulation is 
exact so far. The evaluation of the series is simplified 
by the Wick's theorem-^ (see Appendix IB")) . Each series 
term in Eq. is of the form X^X±, from Eqs. ©, a 
product of several photon operators A t andA^ , which the 
Wick's theorem resolves into a sum of terms composed of 
a normal product and a number of pairs of contractions. 
The matrix element of each normal product between two 
coherent states is simply products of scalars from the sub- 
stitution of cik — > ctk acting on the coherent ket vector 
and a\ — > a* k acting on the bra vector. The contractions 
from Eq. (fBl>|) are 

(A\A 3 ) = 0, (10) 
(AiA]) = K(t t - tj) = |<? fc |V A ^-<A (11) 

Fig. [1] illustrates the diagrammatic representation of 
each term in the scries expansion of the transforma- 
tion matrix and some partial summations of subseries. 
Fig. QJa) shows a typical term. The rules are: (i) the 
initial and final TLS states represented at times and t, 
(ii) the interactions with the TLS denoted by dots in the 
counterclockwise loop for the time-ordered photon oper- 
ators from time — > t and then in anti-time order back 
t — > 0, (iii) the appropriate TLS states between dots 
labeled as ± and (iv) all possible contractions (dashed 
lines) either on the same time line or between the op- 
posite time lines. Fig. [ljb) shows the only two possible 
types of contractions because of Eqs. (|10lll[) . If the two 
ends of a solid line segment have the same (opposite) TLS 
states, the segment is dressed by an even (odd) number of 
photons whose series sum is depicted by a double (triple) 
line, see Fig. HJc). Fig. QTd) provides an example of the 
transformation matrix with a single contraction 

and dressed states. We stress that these diagrams are 
first order in contraction, but infinite order in the coher- 
ent interaction. 



C. The key results 

The quantum treatment of the time evolution of the 
TLS in terms of a series of photon contractions yields 
a number of notable results. The series expansion in 
the TLS-photon interaction permits identification of the 
physical processes. Well-known are the all-dots terms 
uninterrupted by dotted lines, as in the first term of 
Eq. (|B1[) whose sum yields the Rabi rotation and the 
complete contraction pairs in the last term of the same 
equation whose sum yields the relaxation due to the elec- 
tromagnetic vacuum (see Appendix [C]). In the mixed 
terms, the dots between the ends of contractions may 
still be summed as coherent dynamics of the TLS. The 
summation is necessary because along the dressed line, 
each photon gives a term of the order c/|ev|io for a pulse 
of duration to- These series expressions can be summed 





FIG. 1. Diagrammatic representations of perturbation terms 
of the transformation matrix p SJ ,, s ' ;SjS '(t, {&})■ (a) An ex- 
ample of six photon operators with two contractions. The 
two short lines labeled 0, t are the time limits of the inte- 
grals. Each dot labeled with a time variable t n represents 
one photon operator, the lower solid arc being time ordered 
and the upper arc anti-time-ordered. A dashed line stands for 
one contraction between two photonic operators. The uncon- 
tracted dots form a normal product of the photon operators 
for the matrix element of the coherent states. The ± sign 
denotes the state of the TLS at different times, (b) The only 
possible contractions drawn between photons on the same or 
opposite time lines, (c) The dressed line by a sum of all even 
or all odd numbers of photons interacting with TLS without 
contraction, (d) Three diagrams that contains only a single 
contraction for the transformation matrix p++ ; ++. They are 
the leading contributions to the control noise. 



exactly and are given in Appendix [C] The coexistence 
of the coherent dynamics and relaxation yields an effect 
which is extra to the sum of the two processes, as will 
be shown next. This is also clear in the nature of the 
contraction between the evolution and its inverse shown 
in Fig.QJa) and (d). Note how, without the manual sep- 
aration of the photon Hamiltonian into a control part 
and a bath part, the method produces the dissipation 
and the Rabi rotation. More importantly, it contains 
quite explicitly the interference effects between the two 
processes. 

A useful result for our stated purpose of studying the 
fidelity of the control process to high accuracy is the find- 
ing of a small parameter for the expansion. After the con- 
struction of the dressed lines of the coherent processes, 
the expansion of the contraction functions is a perturba- 
tion series in powers of (£o/2i) 7 in terms of the operation 
time to and the decay time T\ (see Appendix[C]for the de- 
coherence time T2 = 2 1 / 7 Ti). The decay in each contrac- 



4 



tion line is ~ / / dUdtjK(U ~ tj) ~ O[(i /Ii) 7 ], where 
the parameter 7 between (1,2) depends on the shape 
of the photon DOS. Therefore, this method provides an 
excellent evaluation of the control noise in the regime 
to <c Ti, while the area of the pulse ~ <?|a|to ~ 0(1). 
The one-contraction terms in Fig. [Tf d) together give the 
consistent result to the first order contribution of the con- 
trol noise which will be evaluated next. The consistency 
of the three diagrams comes from the differentiation of 
a self-energy diagram with one contraction resulting in 
three terms, as in the Ward identity in field theory^ 
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III. RELAXATION AND CONTROL FIDELITY 

The quantum effects between the coherent control and 
the vacuum decoherence depend crucially on the contrac- 
tion, which, Eq. (fTTj) shows, depends on the materials 
properties of the TLS and the photon confinement in the 
form of the DOS of the photonic field weighed by the 
interaction mode dependence. To show explicitly the de- 
pendence of the control dynamics on the weighted DOS, 
it is convenient to model it as a Gaussian: 



p(w) 



k 



g 2 T c 

2^' 



f(w-wo) /4 



(12) 



where g is the average coupling strength and r c describes 
the correlation time of multimodc light in the presence 
of a TLS. The qualitative results, such as the concept 
of the exponent 7 in decay, are unchanged for a general 
DOS that includes a correct behavior for u> — > + . When 
r c — > 0, the broadband DOS yields T\ = l/y/irg 2 T c ; while 
as t c — > 00, the single mode scenario pertains. Under 
this Gaussian DOS, Eq. (|C1|) shows that the contraction 
function is also a Gaussian, i.e., K (t) = g 2 e~ l ■ 

For a broadband DOS, e.g., in free space, the system 
is Markovian and the result is equivalent to that from 
solving the optical Bloch equations.— For an extremely 
narrow DOS, single-mode cavity quantum electrodynam- 
ics dominates. In this section, we use a Gaussian DOS 
whose variable width causes the decay of the upper state 
of the TLS to have a dependence of In P++(i) oc — (i/Ti) 7 
characterized by the exponent 7 and use the change of 
the system from the exponential decay (7 = 1) in the 
broad DOS limit to a Gaussian decay (7 = 2) in the 
narrow DOS region to show the emergence of the quan- 
tum effects of interference between the laser control and 
the vacuum decoherence. The interference actually has a 
beneficial effect on the fidelity of the quantum operation 
on the TLS. 

An example of a relevant physical system is a multi- 
mode finite Q cavity system with bandwidth ~ 1/io- For 
a photonic crystal nanocavity^ with Q ~ 10 4 -10 5 and 
uio ~ 10 15 Hz, the estimated value of r c is ~ 100 ps and 
non-Markovian features could be observable for t ~ r c . 
Then, the small parameter condition for a fidelity of 
0.9999 could be relaxed from £ /7i ~ 10~ 4 to 10" 2 . 



FIG. 2. (color online). Relaxation of an initially excited TLS 
state in vacuum using a Gaussian photon DOS with gr c — 0.1. 
The decoherence of the TLS (solid curve) is different from the 
purely Gaussian (dashed) or Markovian (dash-dotted) decay. 



A. Non-markovian relaxation for a Gaussian DOS 

The short time decoherence behavior of the non- 
Markovian multimode JC system in vacuum is given by 
diagrams with only one contraction (see Fig. [7f b)). For 
an initially excited TLS, we have 



\n[P ++ (t)} ps - 2 / dh dt 2 K{h - ta) 



(13) 



An illustration using gr c = 0.1 is presented in Fig. [2] 
The result reveals that for a finite value of r c . the sys- 
tem evolves from a non-Markovian Gaussian dependence 
(well-known in the onset of the collapse and revival 
phenomena , 14 i 17 ~ — ), P ++ (t) ps , to the Markovian 

exponential decay, P++(t) ps e~*/ Tl . The crossover takes 
place at t ~ r c . For short time, the contraction function 
is almost flat and the decay resembles the single mode 
case. As time increases and exceeds r c , the Gaussian con- 
traction function approaches the broadband Markovian 
limit. Note that both the Markovian and single mode ap- 
proximations overestimate the decoherence of the TLS. 

Fig. [3] depicts the evolution of P ++ (t) in the log scale 
for different gr c . It shows the same crossover from the 
non-Markovian to Markovian relaxation when time is 
comparable to the correlation time. The decay exponent 
7, defined by P++(f) ps e _ ^/ Tl ^ T can vary from 2 to 1 as 
time increases. Because of this non-Markovian dynam- 
ics, the fidelity of the TLS under a coherent control also 
shares the same feature. We will investigate this more 
quantitatively in the next subsection. 



B. Control fidelity analysis 

The control noise problem is solved by the diagram- 
matic method in Sec. [II] with TLS initially in the ex- 
cited state and under a coherent 2tt pulse with truncated 
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FIG. 3. (color online). Plots of log 10 [ — In P++(t)] in short 
time for different values of correlation r c . Crossover between 
Markovian and non-Markovian relaxation occurs at t ~ r c . 
The inset shows the decay exponent, 7. 



Gaussian shape denned by f2(t) = f2 e _( - 2 ' /<T for 

< t < to and zero otherwise. The coherent state ampli- 
tudes a k is related to Vt(t) by Eq. (fC3|) . We take cr = t /4 
and use the Gaussian photon DOS given by Eq. (fT2|) . 

We evaluate the single contraction diagrams in 
Fig. QJd) for decoherence of ~ 0[{t/Tr) J ). The fidelity 
is computed by F(t = t a ) = Tr [P id eaiP(t = t Q )], where 
Pideai is the ideal reduced density matrix. The error, 

1 — F(t = to) is plotted as a function of control dura- 
tion gto in Fig. |4] and is compared with the Markovian 
(i.e. broadband) approximation. The third diagram in 
Fig. |TJd) gives a contribution from the contraction be- 
tween two time lines which tends to be opposite in sign to 
the dissipation effect of the two graphs with contractions 
within the the same line. This is evidence of quantum in- 
terference between control and dissipation. Owing to the 
non-Markovian relaxation, the error of the operation goes 
quadratically with to for small to and then becomes lin- 
ear in to , which is different from the linear to dependence 
in the Markovian limit. The crossover occurs at t ~ t c 
for the same physical reason in the vacuum relaxation 
process. The fidelity is in the form, F = 1 — c(fo/21) 7 , 
where I < 7 < 2 and c is a constant. 

The relaxation at t > to after the pulse finishes is not 
as simple as the Markovian limit, and is given by 

Paa' (t > t ) = Poo> (t = t ) - faa' [A(t)] X [(t - to)/^] 7 . 

(14) 

The renormalization factor f aa i [Ait)] is a functional of 
A(t), depending on the shape of the pulse. This reflects 
the history-dependent dynamics of the TLS when we con- 
sider a Gaussian photon DOS. Such a feature stands in 
contrast with the Markovian case, where the function 
faa'[A(t = t )] only depends on the total area of the 
pulse. 

The discussion above is based on a pure starting state. 
Consider the extension to a mixed starting state. Let the 
prepared pure state at t = relax to the mixed state at 
time U > when a control operation starts. Then, the 




FIG. 4. (color online). A plot of error, l-F(t = to) for a TLS 
driven by a nominal 2ir Gaussian pulse using a Gaussian DOS 
with gr c = 0.03 (solid blue) and a constant DOS (dashed red). 
For the Gaussian DOS, the to dependence of error changes 
from quadratic to linear when to ~ r c , whereas the error 
depends linearly on to in the Markovian regime. 



TLS is driven by the control pulse and influenced by the 
vacuum decoherence. The theory is the same except the 
laser pulse is shifted in time to O(t) ^ in t j < t < tf. 
The fidelity of an example is studied, using a 2tt pulse, 
which is Gaussian cut off at both ends of the same time 
interval a from the center as before. The results are given 
in Fig.[5]for the relative error defined as the proportionate 
change of state fidelity between the start and the finish 
of the pulse, Rf = 1 — [F(tf)/F(ti)} of the operation as a 
function of the initial time of the pulse. For a small and 
fixed duration to = t/ — U, we observe that the relative 
error increases linearly with t; and then saturates to a 
constant after tj «~ t c . This can be roughly understood 
by employing the decoherence crossover picture discussed 
in subsection If If Al When the operation is performed 
inside the non-Markovian region, F - e -^ 2 and R F ~ 




0.000 0.005 0.010 0.015 0.020 0.025 0.030 

FIG. 5. (color online). Relative error, R F = 1 - F(t f )/F(U) 
as a function of ti of a TLS under a 2ir Gaussian pulse for 
fixed gto — g(tf — ti) = 0.01 and gr c — 0.01. For a Gaussian 
DOS, a change of dependence on ti occurs when the system 
evolves from the non-Markovian to Markovian regime. 
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FIG. 6. (color online) AP++(t), the deviation of P++(t) from 
the exact solution of a driven single mode JC system using 
different methods, n = 1007T 2 , corresponding to a 4-7T rotation 
at gt = 0.2. The errors of the ME approximations are as 
much as that of the classical solution. 



~ 0[(gt) 4 ] and the actual computed results within an 
error of 10 -4 . This error can be further improved to 
0[(gt) 2 ^ n+1 ^] by including higher order diagrams with n 
contractions. The Born approximation and the low or- 
der NZ give results close to the classical treatment of 
the electromagnetic fields can be interpreted as a lack of 
quantum content. The agreement between the low or- 
der TCL and the additive approximation might be an 
indication of the additive nature of TCL. All four ME 
approximations fail to reproduce the quantum interfer- 
ence effects between control and dissipation which are in 
the field theoretic treatment. Note that the problems of 
the NZ and TCL methods are not due to the second or- 
der approximation we used because their errors are also 
of the same order. The failure of the ME methods to 
account for the high accuracy of control fidelity required 
by quantum information processing is illustrated by this 
comparative study. 



2g 2 toti] whereas in the Markovian limit, F ~ e ti ' Tl and 
Rp ~ t /Ti. 

IV. ERROR CHECKS 

We check our method against the exact solution^ of 
the driven single-mode JC model for the non-Markovian 
effects involving control and decay of the TLS and com- 
pare with approximations in the master equation ap- 
proach. It is a limit of the multimode problem for 
gr c = oo and 7 = 2. We choose four commonly 
used ME approximations, (i) the Born series, (ii) the 
Nakajima-Zwanzig (NZ) projection method, (hi) time- 
convolutionless (TCL) projectio n 13 ' 20 , and (iv) the addi- 
tive assumption^ Unlike our error bound, these theories 
contain none. Their results are compared with the exact, 
the diagrammatic and the classical Rabi solutions. For 
details of how the first three methods are used in the cal- 
culations, see Appendix [D] The additive ME, assumes 
that the driving terms of the equation of motion by con- 
trol and by dissipation simply add: 

^p.(t) = -i[H cl (t),p s (t)]+k Ps (t), (15) 
at 

where p s (t) is the reduced density matrix of the TLS, 
H c i(t) the classical control Hamiltonian, and K an ap- 
propriate super-operator on p s (t) for the population de- 
cay and decoherence effect in the absence of a control 
(see Appendix [D|) . 

The TLS is initially in the excited state and driven by 
a single mode coherent state at resonance with the TLS. 
The average number of photons, n = 100 7r 2 , is chosen 
so that, at gt = 0.2, the area of the classical pulse is Air. 
Fig. M plots the deviations of the upper state population 
from the exact solution, AP++(t) = P++(t) - P^ ct (t), 
for the six methods above. The diagrammatic method is 
indistinguishable from the exact solution in the small pa- 
rameter region of gt < 0.2 with a theoretical error bound 



V. CONCLUSION 

We have developed a rigorous solution of the multi- 
mode Jaynes-Cummings model for a two-level system un- 
der the control of a coherent electromagnetic pulse. It 
treats the control process exactly within a rigorous error 
bound for the controlled dynamics of the entire system 
of the small quantum object and the photons to any or- 
der in the small parameter in operation time versus the 
decoherence time, (io/^2) 7 where 7 varies from 1 to 2 
as the density of states of the photons varies from flat 
as in free space to a sharp peak as in high Q cavity. 
Our theory is quantum in the sense of no stochastic as- 
sumption and treats the entire quantum system correctly 
within any given error bound. A diagrammatic represen- 
tation provides a simple picture of different physical pro- 
cesses, including the familiar limiting cases of coherent 
Rabi oscillation and vacuum decoherence. It pinpoints 
the control-dissipation interference as the quantum ef- 
fects found by works which precede ours. The effects are 
relevant to basic quantum phenomena and to technolog- 
ical applications. 

This work shows that vacuum relaxation is not the only 
physical process that results in quantum noise in non- 
Markovian systems. Fundamental quantum noise due to 
interference between the control field and relaxation ex- 
ists and it leads to a decoherence comparable to the vac- 
uum fluctuations. Moreover, our method is not restricted 
to the Markovian limit and is valid for an arbitrary DOS. 
Thus, it necessarily goes beyond the optical Bloch anal- 
ysis. 

The time evolution study shows that the relaxation of 
the TLS can vary from a Gaussian decay to an exponen- 
tial one, depending on the ratio of time to the correlation 
time that characterizes the photon DOS. In consequence, 
the fidelity of a general single qubit operation on a pure 
state can be cast in the form F = 1 — c^o/T^) 7 for a con- 
stant c. The Markovian approximation (7 = 1) ovcresti- 
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mates the decoherence of the TLS. The Gaussian deco- 
herence (7 = 2) provides a lower bound for the quantum 
error of a light controlled TLS. 

The field theoretic technique provides a completely 
quantum mechanical description of a small quantum sys- 
tem interacting with the photonic environment. The con- 
trol noise issue, being important in quantum computing, 
serves as an concrete example to demonstrate its capabil- 
ity. This approach may also be applied to other systems 
under nonclassical photon states. Besides the quantum 
object, the field theoretic technique also permits a calcu- 
lation of the physical quantities of the environment, e.g. 
the quantum feedback on the electric field, photon cor- 
relation functions, correlation between the TLS and the 
photons, etc. 

Comparative studies with the existing master equation 
approximations show their general lack of the quantum 
effects due to the control-dissipation interplay, which is 
not restricted to the single-mode model tested. We hope 
that the results of our rather cursory study of these ap- 
proximations would stimulate more developments in rig- 
orously bounded approximations for the master equation 
and further understanding of quantum effects by the con- 
trast and complementarity between the master equation 
and the field-theoretic approaches. 

Our current theory only considers a coherent state with 
a constant phase. In future work, this can be generalized 
to describe an ensemble of coherent states with a mix- 
ture of phases in order to evaluate the phase error as an 
extension of previous single mode study^ 



ACKNOWLEDGMENTS 

This research was supported by the U.S. Army Re- 
search Office MURI award W911NF0910406. We would 
like to thank Paul R. Berman and Renbao Liu for helpful 
discussions. C. K. Chan thanks Wen Yang for a useful 
comment. 



VnwA({g'k}) i s the TLS-field coupling in RWA with gk re- 
placed by g' k = 2gfcWo/(wo + uJk)- Therefore, the counter- 
rotating terms leads to an effective perturbation V' as a 
correction to the RWA. It is an order O(ga/uio) smaller 
than Vrwa- Hence, in the regime where l/u>o <C ~ 
t <C T\ , the effect of the counter-rotating terms is negli- 
gible compared with decoherence. 



Appendix B: Application of the Wick's theorem 

The evaluation of the general perturbation term is sim- 
plified by the Wick's theorem, a general bosonic operator 
defined byW = Y[ i O i , where L = Y,k( u kih + Vkiby is 
a linear combination of bosonic operators bk and b k , and 
can be rearranged as*i£ 

W = :W:+J2 -- W n,n ■ (OiM 

+ E ■■Wi 1 , jM :(O il O jl }(O i2 O h ) 

O'l ,Jl)#(«2,j2) 



E (O il O jl )....(O in O j 



(Bl) 



(il,2l)...(»n,Jn) 



where : W : is the normal ordered form of W, defined 
by all the creation operators to the left of the annihila- 
tion operators; :Wj,k,l,...' is the normal ordered form of 
Yii^j k 1 Q m wn i cn the operators Oj etc. are left out; 
and the contraction between Oi and Oj is defined only 
for i < j by a scalar, 



{(),(), = <>,(), ■(),() 



(B2) 



Note that the suffix in the O operator denotes its order 
in the W expression, rather than the time index in the 
photon operator A\. 



Appendix A: Rotating wave approximation 

To prove its validity for the control noise problem, we 
start with the Hamiltonian H = Ho + V in Eq. (JTJ with 

the full interaction term, V = J2k 9k{o\. + ak)<J x - To this, 
we apply a unitary transformation e , where^ 



*=E 



9k 



UJ k + U 



{a\ + a fc )cr 2 



(Al) 



Appendix C: Evaluation of control and dissipation 



This supplements Sec. lIIII with an analysis of the math- 
ematical structure of the diagrammatic series, two exam- 
ples of exactly soluble situations, and more details of the 
evaluation of formulas. 



1. The contraction formula 



For g\a\/uio <C 1, we expand the transformed Hamilto- 
nian H = e s He~ s up to second order in gk as: 



H = H Q + V RWA ({g' k }) + V, 



where V' 



7— Cz l~ 



EsU4 

k 



a-k. 



(A2) 
(A3) 



Central to the vacuum decoherence is the contraction 
function from Eq. pip given by, 

{A{t)A{t')^) = K(t -*')=/ rfwp(w)e i(w °- w)(t -''l;pi) 

where t and t' can be on the same or opposite time lines. 
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(a) 






+ + iV, r<7 + 



' \i w - 





(C) 



FIG. 7. (a) gives the diagrams for the classical Rabi solutions 
for p++ i ++(t) and p — (b) and (c) provide all the dia- 
grams in the absence of control for p++ t ++(t) and p_+._ + (t) 
respectively. In the Markovian limit, the two ends of a dashed 
line are squeezed to the same point, (b) consists of two low- 
est order diagrams while there is only one in (c), implying 
Ti — T2/2 in the Markovian regime. 



2. Dressed TLS state lines 

A time segment between any two consecutive points 
chosen from 0, t or contracted interaction points in a di- 
agram in Fig. [1] may be viewed as a dressed TLS line in 
which all possible interaction terms for the coherent state 
matrix element are summed. An example term is given 
by: 



it 



*2re+j + l 



dtzn 



+J 



dtj+2 / dtj + i 



x ({a}\ : A 2n+j Al n+j _ v ..A j+2 At +1 



\{a}) (C2) 



The normal product in the integrand enables the eval- 
uation of the coherent state matrix element by putting 
the creation operator to the left with the substitution 
at — > a* k , and the annihilation operator to the right with 
a/,. — > afe. The control is given by, 

k 



e - l ^° t+ rt +c.c.,(C3) 



which relates ctfc with the envelope function fi(t), the res- 
onant frequency uj and the phase (f> of the driving field. 
Thus, an infinite series sum of integrals like Eq. (|C2[) can 
be carried out, leading to even and odd types of dressed 
photon lines (represented by Fig.QTc)): 

D e (t, t>) = cos Q(t _ 0> 



D 



a (t, t 1 ) = (±)i sin 2 A{t>) ^ 0(t - t') (e**), 

(C4) 



corresponding, respectively, to the double and triple lines 
that are dressed by an even and odd number of photons 
in the control. Note that the dressed function D (t,t') 
picks up a + (— ) sign, when the triple line is on the up- 
per (lower) time line, and gains a phase e 1 ^ (e~^) : if 
the triple line goes from — to + in the clockwise (an- 
ticlockwise) sense. Here, A(t) = /„ dt'fi(t') gives the 
area of the pulse at time t. In the absence of control, 
D e {t,t') -> e(t-t') and D (t,t') -> 0, so that the double 
line is reduced to a single line and diagrams that contain 
any triple line vanish. 

The mathematical expressions representing different 
diagrams can be obtained by first writing down all the 
contraction and dressed functions, and then integrating 
over all time variables corresponding to the vertices of 
each contraction line. The vertex picks up a factor of 
i (— i) if it is on the upper (lower) time line. In the fol- 
lowing, we will provide some explicit examples. 

The classical Rabi solution corresponds to dia- 
grams with no contraction. The two dressed lines in 
Fig. [TJa) gives the transformation matrices p ++ . ++ (t) 
and p !_;++(£). Using Eqs. (|C4[) . we have 



A(t) 



_ ie i(uj t+4>) cog 



A(t) 



A(t) 



.(C5) 



which is the Rabi solution without decoherence. 

3. Vacuum relaxation of TLS 

In the absence of control, the photon lines are not 
dressed. A broadband DOS yields the Markovian limit, 



1 



K M {t-t') = -8{t-t') 



(C6) 



Then all the diagrams in Fig.fTJb) and (c) can be summed 
exactly, yielding respectively: 



vacuum, M / ,\ 

p++,++ (*) 



vacuum. M / ,\ 

P-+.-+ (t) 



o iu t-t/2T 1 



-+,-+ w ~ - (C7) 
This Markovian limit leads to the standard result of spon- 
taneous emission, where T 2 = 2T\. This relation can also 
be seen from the lowest order terms in that p^?^ 7 ?' (t) 

vacuum. M 



it) 



contains two lowest order diagrams, while p 1 |_ 

only one. 

Using the same argument, for a non-Markovian sys- 
tem with a general DOS and the 7 parameter defined in 
Sec. IIII Al we have 



vacuum (i\ -i 

P++. + + \ l ) — 1 



vacuum /jA 

P-+,-+ w ~~ 



so that T 2 = 2 1 /T i . 



1 2 



O 



2T 



+0 



27' 
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FIG. 8. (color online). P++(t) « P} u j(t) + P£l(t) of an 
initially excited TLS under a Gaussian pulse in the Markovian 
limit. It is in quantitative agreement with the optical Bloch 
solution. 



4. Exact solution to first order in contraction 

Fig. [Ud) s h° ws a U the diagrams with one contraction 
with the dressed lines. By the diagrammatic rules, the 
transformation matrix is: 

p<& t++ {t)=2Z 1 (t) + Z 2 (t), (C9) 

where Z x {t) = -D e (t,0) [ dt f dt'K(t-t') 
Jo Jo 

x D e (t,t)D e (t,t')D e {t',0), (CIO) 

Z 2 (t)= [ dt [ dt'K(t-t') 
Jo Jo 

x D o (t,t)D o (t,t')D e (t,0)D e (i',0). (Cll) 

Z\ and Z 2 correspond to the non-crossing and crossing 
diagrams respectively. This result is valid for arbitrary 
photon DOS and thus covers both the Markovian and 
non-Mar kovian regimes. This is the basis for the result 
in Sec [TITI using a Gaussian DOS, Eq. ([T2"]). 

The diagrammatic method can reproduce the Marko- 
vian results. By Eq. (|C6)) and ([05]) . Fig. U shows a 
Markovian example of an initially excited TLS driven 
by a Gaussian pulse at resonance (the same physical sit- 
uation as Sec. IIIII B). Our result not only shows an ex- 
cellent agreement with the optical Bloch analysis in this 
Markovian limit, but also allows an understanding of the 
underlying processes. For instance, the nodes of P ++ (t) 
do not reach zero. We note that this feature has been 
observed in many experiments but cannot exclude 
other experimental noise sources as the cause of the fea- 
ture. This theoretical feature can be understood through 
our diagrammatic representation in Fig. [TJd) . One can 
show that the non-zero node arises from the crossing con- 
traction term, which is very different from the vacuum 
decoherence that contains non-crossing contractions only 

(Fig. ED)))- 



Appendix D: Master Equation Approximations 

The four approximations for the ME approach used in 
Sec. HV] are detailed here. The standard ME for the TLS 
up to the second Born approximation is given by (p. 250 
of Ref.IH): 

^ Born 1 1 \ 

dt p ° {t) 

= -iTr R [V(t),pf°™(0)®p R ] 

dt'Tr R [V(t), [V(t'),p?° rn {t') ® PR \] . (Dl) 

This equation is derived from a second order expansion 
of the Liouville equation. In the single mode JC system, 
at resonance, the interaction is V(t) = g (er+a + er_a^) 
and the reservoir density matrix constant, p R = \a) (a\. 

The Nakajima-Zwanzig (NZ) and time convolution- 
less (TCL) projective operator techniques are outlined 
systematically 13 i 20 The usual assumption of Ti^ [V(t) <g) 
pa] vanishing is unnecessary and is not made in our con- 
trol problem. The second order NZ ME, derived by the 
method of Brcuer et a/. i 13 ' 20 is then: 



dt 



iTv R [V(t),p* z (t)®p R ] 

dt'Tr R [V(t), Tv R [V(t'), p™(t') ® p R ] ® p R ] 

dt'Tm [V(t), [V{t% p^ z {t') p R }] . (D2) 



The second order TCL ME shares the same structure but 
p^ z {t') is replaced by p^ GL (t) in the integrand. 

The additive ME assumes that the control and relax- 
ation terms are additive (see Eq. (fTS")) ). In the single 
mode JC system, it becomes: 



(D3) 



,-i<t> 



= - ig\a\[a + e + a_e i0 , p| dd (t)] +gta,ngt 

x {2a„pf d (t)* + - a+a-pf d {t) - pf d (t)a + ^} , 

where a = |a|e 8 ^. By neglecting the second term, the 
first control term leads to the classical Rabi motion. On 
the other hand, in the absence of the control (\a\ = 0), 
the second term will produce the vacuum Rabi oscilla- 
tions of the single mode JC system. In the small time 
regime (gt <^ 1), it corresponds to a Gaussian relaxation. 

For an initially excited single mode JC system under a 
coherent control with A = 2g\a\t - O(l), Eqs. (|DHID3p 
are solved and compared with the exact, the diagram- 
matic and the classical methods. The diagrammatic solu- 
tion is calculated from Appendix [Cl using the single-mode 
contraction function K(t—t') = g . The comparison uses 
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the resultant solutions: 

Irvl 2 ™ 

pexact ( ^ = ^cos 2 (gV^TT t) 

n 

P^ gram (t) = cos 2 5 |a|i - M. cos2.g|a|t 
3 fft 



■ sm . 

lal 



i2g\a\t + 0[(gt) 4 ] , 

pclassical^ = ^2 ^| a |^ _ QyQ 
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